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Abstract. We investigated how the presence of an additional lattice potential, driven 
by a harmonic noise process, changes the transition rate from the ground band to the 
first excited band in a Wannier-Stark system. Alongside numerical simulations, we 
present two models that capture the essential features of the dynamics. The first 
model uses a noise-driven Landau-Zener approximation and describes the short time 
i evolution of the full system very well. The second model assumes that the noise 

CNl ■ process' correlation time is much larger than the internal timescale of the system, 

yet it allows for good estimates of the observed transition rates and gives a simple 
interpretation of the dynamics. One of the central results is that we obtain a way to 
control the interband transitions with the help of the second lattice. This could readily 
\ be realized in state-of-the-art experiments using either Bose-Einstein condensates or 

CO ' optical pulses in engineered potentials. 



1. Introduction 



During the last years, there have been enormous advances in the experimental control 
over Bose-Einstein condensates (BECs) in optical lattices. States can be prepared with 
a high degree of control [T] and the systems can be isolated from external perturbations 
for long timescales [2]. This fine control over the system allows to investigate quantum 
mechanical effects that were previously very difficult to observe in experiments. Among 
these effects are Bloch oscillations, Wannier-Stark ladders, and Landau-Zener tunneling 
[3H5]. Dynamics similar to the ones observed in dilute BECs can nowadays also be 
realized using purely optical means [HIE]- 

These advances have also paved the way to studying transport properties of spatially 
[8HT2] or temporally disordered systems [13J. Disorder is to some degree present in 
almost all naturally occurring systems; in particular it plays a large role in solid state 
systems. New effects caused by disorder, such as Anderson localization have therefore 
potential ramifications in a wide range of fields. 
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But systems influenced by noise are not only an interesting subject due to their 
omnipresence; a noise process can also allow to fine-tune the properties of the system, as, 
for example, observed in the effect of stochastic resonance |14j . The question naturally 
arises, whether a similar effect can be seen in an optical lattice system driven by a noise 
process. 

While most investigations of noise-driven systems focus on white [15] or 
exponentially driven noise [T&H20] . we use a process known as harmonic noise in this 
work. This noise process exhibits a strongly peaked power spectrum, corresponding 
to an oscillatory behaviour. It therefore introduces a characteristic noise frequency 
and has been observed to facilitate tunneling processes if this frequency matches the 
characteristic frequency of the system [2T1 [22] . In stochastic dynamics, it is now 
understood that the escape between attractors takes place along well defined "paths" in 
phase space, which are typically analytic [T8~ tl2"3H2"6] : it is the matching between these 
structures and the stochastic noise which facilitates the escape, a phenomenon observed 
also in transient dynamics [27]. 

The system we consider in this work consists of a static tilted optical lattice which 
is disturbed by a second, noise-driven lattice. In the absence of the noisy lattice, the 
system corresponds to the well-known Wannier-Stark problem [28]. This Wannier-Stark 
system will be used as a reference system throughout this work. 

Previous work by Tayebirad et al. [29J reports results based on numerical 
simulations of our system. In this paper, we will investigate the effects of the parameters 
of the noise process more thoroughly and introduce two different models that indeed 
explain many facets of the behaviour of the noise-driven problem. 

In section |2] we will introduce the system studied in this paper and discuss it in the 
context of non-interacting BECs in optical lattices. We then introduce two models that 
allow for an interpretation of the behaviour of the system in section |3j An extensive 
numerical study is then presented in section HJ these results are compared to the models 
presented in the previous section and we give an outlook on weakly interacting BECs. 
We close with a short conclusion of the results in section |5j 



2. Noisy Wannier-Stark system 

In dimensionless units the Hamiltonian studied in this work reads 

H = - l -d 2 x + V [cos(x) + cos(a[x - 0(t)])] - F x, (1) 

where Vn gives the potential depth of the lattice, F is the Stark force, a is a real 
numbeiyj and 4>{t) is a so-called harmonic noise process (introduced below). For a 
conversion between our dimensionless units, Si-units and common experimental units 



(applied e.g. in [5l[29j|30]), see appendix |Appendix B 



| If this number a is irrational, the resulting system is quasi-periodic, similar to those used to study 
Anderson localization by Roati et al. [12] , 
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In the absence of the second phase-shifted stochastic lattice, the Hamiltonian in 
(PQ) corresponds to the well known Wannier-Stark system [311 [32]. A quantum state 
|^), which is sufficiently localized in momentum space (i.e. much smaller than the 
Brillouin zone) and prepared in the ground band, undergoes Bloch oscillations with the 
period Tb = 1 /Fq and partially tunnels into the first excited band everytime it reaches 
the edge of the Brillouin zone (corresponding to an avoided crossing of the respective 
energy bands) . The characteristic time scale and frequency of the system are thus the 
Bloch period Tb and Bloch frequency ojb = 27r /Tb [31]. The second potential term 
renders the system stochastic (via the phase 0) and spatially disordered in the static 
case (<f) = const) if a ^ Q. 

Possible experimental realizations of such systems are for example BECs in optical 
lattices [HIE]- Similar Hamiltonians are for example used for precision measurements 
of forces acting on the atoms of the BEC [33]. If the BEC is either dilute enough 
and/or interactions are suppressed by tuning with the help of a Feshbach resonance, 
their behavior can be well described by a single particle Hamiltonian as the one in 
(|T|) (see, e.g., [5l[30l[3l])- In order to better understand how the observed effects would 
manifest in BECs with non- vanishing interactions, we will also look at simulations using 
the Gross-Pitaevskii equation in section 14.21 There are also recent experiments using 
optical pulses in engineered potentials that realize similar Hamiltonians [6j[7]. 

Let us briefly summarize the most important properties of harmonic noise. 
Harmonic noise can be imagined as a damped harmonic oscillator driven by gaussian 
white noise. It is defined via the two coupled stochastic differential equations [2T1I22] 

= H (2) 
ft = -2I> - w o 2 + V4FTZ(t), (3) 

where T is the damping coefficient, uq sets the characteristic frequency of the noise 
and T determines the noise strength. represents gaussian white noise [351I3S]- The 
equilibrium distributions of the noise variables are of bivariate gaussian type with first 
and second moments given by [221135] 

((/>) = 0, M = o, (4> 2 ) = ^ and (/^> = T - ( 4 ) 

Just as the damped harmonic oscillator, harmonic noise can be classified into two 
different regimes, Uq > \/2T and loq < y/2T. In this paper we will only consider the first 
case, where the noise shows oscillatory behaviour and its power spectrum has a clear 
peak at a frequency oj\ = a/ w o — ^r 2 « oo for uj ^> T [22 ] [29 ] 135]. 

While the work of Tayebirad et al. provided a first numerical investigation of this 
system and showed that the noise properties have an influence on the observed transport, 
it could not answer the question after the mechanism of this influence [29]. In order 
to provide an answer to this problem, we present two simple models that explain how 
the noise process acts on the wavefunction in different regimes. We then compare the 
predictions of the two models with a more systematic numerical study and find that 
they give a good agreement with, and provide an intuitive understanding of the relevant 
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dynamics. Finally, we also look at the stability of the observed effects against nonlinear 
interactions as they would, e.g., be present in a realization using BECs. 



3. Simplified Models 



3.1. 'Noisy' Landau- Zener model 

It turned out that the incommensurability of the two lattices does not have a great 
effect on the survival probability [37] due to the time-dependent phase (pit). We will 
thus set a = 1 and present a two state model that allows us to approximate the system's 
dynamics around the band edge. The Hamiltonian of ([I]) now reads, 

1 



H 



-d x + V [cos(x) + cos(x — 4>(t))] — F x. 



(5) 



The translational invariance (at fixed time) of this Hamiltonian can be recovered by 
applying the gauge transformation ip = e lFoxt ip (changing the frame of reference from 
the lab system to the accelerated lattice frame) [5j[3T]- This yields, 



H = - (p + F tf + V [cos(x) + cos(x - <f>(t))] 



(6) 



where we identified — id x as the momentum operator p. It is instructive to express this 
Hamiltonian in the momentum basis, this gives 



H 



dp- [{p + Fotf \p)(p\ + V ((1 + e**)\p)(l+p\ + (1 + e"*)|l + p){p\)] , (7) 



from which it is clear that it only couples states with Ap = p — p' G Z (due to 
the time evolution). The momentum states can thus be written as \p) = \n + k), 
with n G Zo and k G [—0.5,0.5) C K being the quasimomentum of the system. Due 
to the conservation of quasimomentum, the Hamiltonian in (J7J can be decomposed 
into independent Hamiltonians each acting on a subsystem of constant quasimomentum 
k [3[38]- A system initially in state \po) = \k + n) thus evolves according to the 
(tridiagonal) matrix Hamiltonian, 



/ 



H, 



ho 



Vb(l + e^) 
(k - 1 + F t) 2 
V (l + e-^) 



0\ 



v 







Vo(l + e^) 
(k + F t) 2 

Vb(l + e"^) 



Vb(l + e*) 
{ko + l + Fot) 2 

V (l + e-'*) 



(8) 



/ 



where the Vo(l + e terms stem from the optical lattices and enters as a complex 
phase. The dynamics of the full system around an avoided crossing of the ground and 
first excited energy band (at the edge of the Brillouin zone) can be approximated by 
the highlighted part of the above matrix^ with k = 0.5, i.e. the value at the edge of 

§ This corresponds to reducing the system to the two lowest energy states. At the edge of the Brillouin 
zone, this approximation is very accurate as long as the off-diagonal coupling terms are not too large, 
since only the lowest energy states contribute to the first excited as well as the ground band. 
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the Brillouin zone jSJEH] • This yields a reduced Hamiltonian 

fr - l ( ~ Fot W + ei M (o) 



where we additionally subtracted \ + (-Fo^) 2 , i-e. we shifted the absolute value of the 
energy scale. This 2-by-2 matrix can be seen as an effective 'noisy' Landau-Zener (LZ) 
Hamiltonian; its instantaneous eigenvalues are given by, 

E N .i,2 = t\ ^(F ty + 2V Q 2 (cos(<P(t)) + l). (10) 

We can thus introduce an effective band gap as a function of time by averaging over the 
probability distribution of the noise process <ft. At t — 0, this leads to 

AE cS = V [(v/2(cos(0(O)) + l)> ± Std (y2(cos(0(O)) + i)))] , (11) 

where (f(4>)) denotes the average over the noise process and 'Std' means standard 
deviation (see figure |6] for a schematic representation of AE e g). 

In the case of a single (noiseless) optical lattice only the constant term in the 
anti-diagonal of the Hamiltonian in is present and one recovers the standard LZ 
model [39H12] . Starting in the ground state at time t = — oo, the probability to remain 
in that state until time t = oo is given by the Landau-Zener formula [39ti42] . 



P sur (t = oo) = 1 - exp ( -^f ) . (12) 



2F 

This allows to give estimates for the survival probability in the 'noisy' LZ model by 
replacing Vq with AE% S in (JT2J). 

The two state approximation around an avoided crossing is only applicable if the 
initial momentum distribution is sufficiently localized in momentum space (i.e. much 
smaller than the Brillouin zone) [30] and if the characteristic time of the system is larger 
than the transition time in the two state model [13JH1], i.e. the transition time should 
be smaller than one Bloch period. 



3.2. Quasistatic model 

Since the noisy LZ model is limited to short timescales and can only be applied if both 
lattices in (pQ) have the same wavelength, we will present another model which is able 
to describe the behaviour of the system in cases where those conditions are not met. 

This model is based on a quasistatic approximation of the noise process; <f>(t) 
is replaced by a term /3t linear in time. The resulting system is analyzed and its 
observables are averaged according to the properties of the noise process. In the following 
paragraphs, we will provide a more detailed description of this model. 

For short timescales, the dynamics of the system can be approximated by a Taylor 
expansion of the noise term (j)(t): 

<f>(t + e) = 4>{t) + e^-Mt) + 0(e 2 ) = <f>{t) + e/i + 0(e 2 ). (13) 
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Figure 1: The fluctuating red curve shows the survival probability for the noiseless 
bichromatic lattice for a relative lattice velocity of 0. The parameters are a = 0.61, F = 
0.00762, T = 0.00762 and V = 0.125. The black bars mark the values of (3 which lead 
to an intertwined position of the barriers (see figure [2]). The blue gaussian curve of 
variance T gives the probability distribution of fi given in (jlj). In the quasistatic model, 
we set = fi and P SVLr {(3) is therefore averaged over this gaussian distribution. 



A rigorous mathematical treatment of the e 2 term in f|T3|) is far from trivial]]]. 
Nevertheless, it is clear that it depends on ^(p(t) = ^l^(t) as given in (J3]) and is 
thus negligible for cuq, T <C 1. In this case, ([I]) can therefore be approximated as 

H = -\d 2 x + V cos(x) + V cos (a(x - cf> - f3t)) + Fx, (14) 

where the noise variable fi(t) has been replaced by a constant term /3. Equation 
( IT4l) describes a tilted bichromatic optical lattice system where the two lattices have 
a constant relative velocity /3. The dynamics of such a system can be seen in figure [TJ 
The fluctuating red curve shows the survival probability of the state in the ground band 
after £ = T# . It can be seen that this survival probability drops sharply for intermediate 
values of/3 (|/3| ~ 0.31 — 1.31, region between the black vertical lines). This corresponds 
to strongly enhanced interband transport for these parameters. 

As in section I3.1[ the key to understanding lies in considering the effect of optical 
lattices on the wavefunction in momentum space. Let us first look at a single tilted 
optical lattice Vocos(aa:). We assume that the timescale of the Landau-Zener transition 
from the ground band to the first excited band is short compared to a Bloch period 
and is thus accurately described by the Landau-Zener model. Due to the static force F, 
a momentum-eigenstate scans through momentum space; once it reaches a momentum 
p = ±|, it undergoes a reflection on the lattice with the probability of P sur given in ( TT2l) . 
This process is schematically displayed in figure |2] (left). For the ground band, the effect 
of the lattice can therefore be represented as a pair of barriers at ±|a (blue barriers in 
left of figure [2]), which trap the state in between them. This momentum range between 

|| As the second time derivative of the noise term <fi(t) is not well-defined, it is necessary to look at a 
value {<t>{t)) which is averaged over a timescale, which is much shorter than the characteristic one of 
the system. 
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the barriers corresponds to the ground band. For the lattice moving with velocity 0, 
the position of these barriers in momentum space is shifted to + 0). Adding up 
the effects of both, static and moving lattice thus results in the presence of two pairs 
of barriers as seen in the right part of figure [2j If the two pairs of barriers are in an 
intertwined position (as shown in the right of figure E]), reflections on the moving lattice 
can help the state escape from the ground band. This can be understood from the 
right part of figure |2j If the state 'hits' the barrier at k = — ^ + /3 from the left (green 
barriers in figure [2]), it can either be transmitted or reflected. While the transmitted 
fraction of the state is still trapped in-between the barriers at k = the reflected 
fraction appears at the right side of the barrier at k = | + j5 and has thus escaped the 
ground band of the non-moving lattice. This can be interpreted as absorbing momentum 
from the moving lattice by means of a Bragg reflection and thus being 'boosted' out 
of the ground band. It should be noted that the same process can also take place for 
negative beta (transporting the state to the left of ground band in figure [2j therefore 
also out of the ground band). Transport to higher energy bands is therefore enhanced 
for intertwined barriers. This is supported by the data visible in figure [U where the 
values of (3 which lead to an intertwined configuration are the area in-between the black 
vertical lines. These values of (3 do indeed correspond to a heavily suppressed survival 
probability in the ground band and thus to strongly enhanced transport. 

Keeping this picture in mind, we can understand how the noise process acts upon 
the momentum states of the system. The constant /3 is responsible for the relative 
position of the barriers in momentum space. In the real system, /3 is given by the noise 
variable fi(t) and is thus a function of the time. The positions of the barriers of the 
noise driven lattice thus move in momentum space. 

In our quasistatic model we assume fi(t) to be a constant /3, which is a realistic 
approximation as long as the noise process is slow, i.e. fi(t) changes slowly with time. 
For n(t) = (3 the survival probability in the ground band is given by figure [TJ To obtain 
an approximation for the behavior of the full, noise-driven system, we thus average this 
survival probability P(/3) over the equilibrium distribution of fi(t), which is a gaussian 
of width T (see (j4|)) as pictured in figured! 

While we would expect this model to be accurate for slow noise processes, it ignores 




Figure 2: Left: Effect of a single lattice on momentum eigenstates. Right: Combined 
effect of a static and moving lattice. 
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the influence the non-static nature of n(t) has on the Landau-Zener transitions and does 
also not account for time-correlations of the noise-process. Especially for a fast noise 
process where fi{t) changes significantly during one Bloch period, we therefore expect 
discrepancies between the quasistatic model and the full system (35]. 

4. Numerical Results 

We compare now the predictions made by the two models introduced in the previous 
section with data from the full system described by the Hamiltonian in ([1]). 

The most important property of the noise process <j) is its variance. In the following, 
we will analyze the survival probability of the system in the ground band for the two 
cases: Var(0) = constant (section H~T]) and Var(0) = changing (section [472]) . The latter 
is realized by keeping uo constant and varying T. 

Figure [3] shows the survival probability in the ground band after one Bloch period. 
The survival probability strongly depends on the noise parameters T and ojq. 
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Figure 3: Survival probability of the 
system described by ([T| in the ground 
band (colour coded) after t = Tb, i.e. 
after one transition, versus the rescaled 
noise frequency ujq/ujb and the noise 
strength T/oj b . The vertical solid line 
represents a line of constant uq. The 
inclined dashed line represents a line of 
constant variance. Data from numerical 
simulations of the Hamiltonian in (151), 



for details see appendix |Appendix A 
The parameters are F = 0.00762, T = 
0.00762 and V = 0.125. 



4-1. Constant variance, i.e. T/u^ = const . 

Here, we show the survival probability in the ground band after one Bloch period for 
constant variance (cut along the dashed line in figure [3]) versus the noise frequency. We 
focus on the region of the initial decrease and minimum of the survival probability, i.e. 
ujq/ujb < 14. Figure HI compares numerical data from the 'noisy' LZ model in (Q, with 
data from the full Hamiltonian in (|5]) for two values of the potential depth Vq. 

There is an excellent qualitative and a good quantitative agreement between the 
'noisy' LZ predictions and the simulations of the full system. The 'noisy' LZ accurately 
reproduces the initial decay of the survival probability (as obtained for the stochastic 
Wannier-Stark system) for increasing noise frequency and gives a good approximation 
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Figure 4: Survival probability in the ground band at t = Tb versus the rescaled noise 
frequency ujq/ujb- Shown are data for the real system (blue dashed lines), the 'noisy' LZ 
model (red solid line) and the quasistatic model (black dotted line). The vertical lines 
indicate the position of the first minimum in the survival probability for their respective 
data sets. The parameters are F = 0.00762, T = 0.00762, (0 2 ) = 0.5 and V = 0.0625 
(left) and V = 0.125 (right). 



for the position of the minimum. Yet, it systematically overestimates the survival 
probability. This overestimation is less pronounced for higher potential depths Vq 
(compare figure EH (left) and figure EH (right)). Most importantly, the position of the 
minimum changes with varying potential depth Vq. The quasi-static model also shows 
a decrease in the survival probability with increasing frequency, but cannot reproduce 
any of the finer features visible in the stochastic Wannier-Stark system and the 'noisy' 
LZ-model. 

The varying position of the minimum leads us to the next figure, figure EJ in 
which the position of the first minimum in the survival probability is plotted versus the 
potential depth Vo. Here, the predictions of the 'noisy' LZ model and the quasistatic 
model are compared to the full system; also the effective band gap AE cfi is plotted 
as a function of Vo- The data for the full system shows that there exists a linear 
relationship between the position of the minimum in 03]) and the potential depth Vq. 
This relationship is accurately reproduced by the 'noisy' LZ model, but it can not be 
seen in the quasistatic model. The fact that the quasistatic model can not account for 
the position of the minimum is easily understood by looking at figure HJ Since the model 
does not predict the local minima visible in P sur of both the full system as well as the 
noisy LZ model, it is clear that the position of the first minimum will in general be at 
much larger values of ujq. 

Another observation that can be made in figure [5] is that the slope of the linear 
curve (for the full system and the LZ model) is approximately given by the one for the 
effective band gap. In fact, linear fits to the data give a slope of 1.93 for the real system 
and 1.76 for the LZ model, compared to 1.88 ± 0.15 for the effective band gap (TTTj) . 
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Hence, within errors, the position of the minimum is given by the effective band gap 
plus a constant off-set. The linear relationship also holds for different values of F [37J. 
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Figure 5: Frequency position of 
minimum versus potential depth 
Vq. Shown are the data for 
the real system (blue lines), the 
'noisy' LZ prediction (red squares), 
the quasistatic model prediction 
(green crosses) and the effective 
band gap AE e g (dashed-dotted line). 
The parameters for the numerical 
simulations are Fq = 0.00762, T = 
0.00762, (<p 2 ) = 0.5. 



The linear dependence of the minimum on Vq can be interpreted as follows. As 
mentioned in section |2j the harmonic noise 'feeds' an energy of ~ loq into the system. 
Once this energy is high enough to overcome the band gap between the ground and first 
excited band, given by ~ A£Lg, the system can be excited into the upper band via a 
'phononic' excitation process^ (red wiggly lines in figure |6]). 

The overestimation of the survival probability by the 'noisy' LZ model occurs 
because the 'LZ' model overestimates the separation of the bands in the real system 
far away from an avoided crossing (the model only works in the blue shaded region 
of figure |6]); hence transitions to higher bands are suppressed (the effective bandgap 
increases linearly with time t for large times, see figure |6]). For large Vo, only transitions 
close to an avoided crossing can occur, but for smaller Vq transitions across the full 
Brillouin zone happen in the real system and this can not be ignored. 




I V^Bloch 



Figure 6: Instantaneous eigenvalues 
(solid black lines) of the 'noisy' LZ 
model for a typical noise realization. 
Schematically represented are AE e g and 
transitions between energy states (red 
wiggly lines). The blue shaded area 
extends over one Bloch period and 
approximates the energy band structure 
of the 'noisy' WSS around an avoided 
crossing. The crossed out transition is 
not allowed. 



% 'Phononic' is not to be taken literally, but indicates that the energy is transferred to the atoms via 
an optical lattice shaking/vibrating with a more or less well-defined frequency of w ujq. 
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The early rise of the survival probability in the 'noisy' LZ model is due to its two 
band nature. If the energy fed into the system by the noise is too large, there exists 
no target state for the transition (see crossed out transition in figure E]). In the real 
Wannier-Stark system, this rise happens only at higher u because it is not just a two 
band system and transitions to higher bands are actually possible. 

4-2. Changing variance, i.e. constant uq and varying T . 

In figure |3] we observe the following: once u assumes a small enough value it stops 
having a discernible influence on the survival probability P sur in the ground band. In 
this section we will therefore take a more detailed look at how P sur in this regime depends 
on T while keeping uq = 1 constant. 

On the left side of figure we see the survival probability of a system initially 
prepared in the ground state in comparison to the predictions from the two models. 
The observable P snr (t = TBioch) is the same as previously looked at in figure HI 
Good qualitative agreement between the full system and both models is observed. 
Furthermore, within the numerical errors, the minimum of the survival probability is 
correctly predicted by both models. There are, however, quantitative discrepancies for 
small (noisy LZ model) as well as large values of T (quasistatic model). 

On the right side of figure [TJ a system with two different lattice constants is studied 
(a = 0.618). The quantity analyzed in this case is the large time survival probability 
of the ground state per Bloch period P sur (t + TB\ oc h) / P S ur(t) for t > 5. This quantity 
is similar to the previously examined survival probability P sur (t = TbiocIi); but it is 




T/ ^Bloch T/ w Bloch 



Figure 7: Survival probability in the ground band for fixed noise parameter ujq. Plotted 
are the results from the full system (blue stars), the predictions by the noisy LZ model 
(red line) as well as the predictions from the quasistatic model (black dashed line). 
Parameters are F = 0.00762, T = 0.00762, V = 0.125 and a = 1 (left) and a = 0.618. 
The plotted quantity is the survival probability after one Bloch period P sur (t = Teioch) 
(left) and the average survival rate for long timescales P SUT (t + TBi oc h)/P S ur(^) (right). 
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derived from the long term exponential decay of the initial state. No results for the 
noisy LZ model are shown since this model cannot make any predictions for the case 
a 7^ 1. For the quasistatic model, good quantitative agreement with the full system is 
observed for most values of T. Nevertheless, significant differences are visible around 
the predicted minimum. These differences, however, result from the limitations of the 
quasistatic model and vanish if uq as well as T are chosen small enough. In this case, 
where q / 1, the quasistatic model is very successful in quantitatively predicting the 
long-term survival probabilities, especially if T as well as uq are small |45j . 
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Figure 8: Survival probability in the 
ground band after one Bloch period 
using the Gross-Pitaevskii equation (same 
parameters as left of figure 0). Results 
from the (linear) Schrodinger equation are 
shown in solid blue, while green crosses 
show results for iV = 5 • 10 4 and red 
stars for iV = 10 5 atoms in a cigar-shaped 
condensate. The simulated setup was the 
one used in the experiments reported in [16] 



and is described in |Appendix A[ The used 



atom numbers correspond to dimensionless 
nonlinearity parameters of C ~ 0.08 (green 
crosses) and C ~ 0.13 (blue stars), c.f. 
ref. \MM- 



In order to understand whether the described effects would be visible in a realistic 
experimental setup using BECs, it is important to know how they are affected by 
interactions between the atoms in a condensate. Using the 3D Gross-Pitaevskii equation, 
we performed numerical calculations for an experimental setup similar to the one used 
in jl6]: 87 Rb atoms that are placed in a cigar-shaped optical dipole trap with 250Hz 
radial and 20Hz longitudinal frequency. The results can be seen in figure |8] and show 
that, while clear differences in the survival probability exist, the enhanced tunneling rate 
for intermediate noise frequencies persists and is hardly influenced by the nonlinearity. 
It should especially be noted that for a nonlinearity parameter (as defined in [46J) up 
to C ~ 0.13, the position of the minimum is hardly changed from the single particle 
dynamics El- 

+ Nonlinearity parameters are only given approximately, as they depend on the maximal density of the 
condensate. Since our setup is governed by a stochastic equation, the peak density of the condensate 
is not only time-dependent due to the interband tunneling [4 7) . but also fluctuates for different noise 
realizations. 
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5. Conclusion 

We have investigated the transport properties of a bichromatic Wannier Stark problem 
under the influence of harmonic noise. Our main interest was understanding the 
circumstances under which the noise process can enhance or suppress the transport 
in momentum space. 

An extensive numerical study of the influence of the noise parameters revealed that 
the system's behaviour can be divided into two different regimes, as seen in figure [3j In 
the regime of low variance of the noise process, the noise only has very little influence, 
while in the regime of large variance, the properties of the system almost exclusively 
depend on the noise strength T. In this regime of large variance, the transport in 
momentum space has a clear maximum for a certain value of T. 

We presented two different models that both offer an interpretation of these results. 
The noisy LZ model approaches the system as a noise-driven Landau-Zener transition 
and argues that transport should be maximized if the noise frequency matches the 
band gap between the ground and the first excited band. While this model can only 
predict results for the case where both lattices have the same wavelength, it reproduces 
the numerical results of the full system in this case very well. Furthermore, it gives 
an accurate prediction of the noise parameters necessary to maximize transport (see 
figure [5]). The quasistatic model assumes the noise process to be slow compared to the 
timescales of the noise-free system. Linking the relative velocity of the optical lattices 
to their action upon the wavefunction in momentum space, this model predicts that 
the transport rate only depends on the noise parameter T and achieves a maximum 
for intermediate values of T. Within the low-variance regime seen in figure [3] this 
prediction matches very well with the results of the full system, even if full quantitative 
agreement is only reached for very slow noise processes. Our quasistatic model, where 
the particle is effectively kicked by the moving second lattice, see figure EJ may turn out 
to be relevant as well for the understanding of damping effects of Bloch oscillations in 
'noisy' solid-state devices [48|l49]. 

Since experiments that study the expansion of the wavefunction in noise-driven 
lattices are underway [13] . an experimental realization of our system and a comparison 
between experimental and numerical data would be interesting. As shown at the end of 
section H~2| the observed effects are fairly stable against interactions between the atoms 
of a BEC and should therefore be easily observable. Besides possible realizations with 
BECs in optical lattices [12] , purely optical techniques such as the one of [6] could also be 
used. While control over the BEC using a noise process is subject to large fluctuations, 
a deterministic bichromatic lattice, as presented in the beginning of section 13.21 could 
be used to control the transport between different bands with a high degree of precision. 
In figured] we, e.g., can see that with small changes of the relative lattice velocity f3, 
the tunneling rates can be easily tuned by various orders of magnitude. 
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Appendix 

Appendix A. Numerical methods 

The data for the 'noisy' LZ model is obtained by numerically integrating the Schrodinger 
equation corresponding to Hamiltonian ((9]). The survival probability has been obtained 
by starting in the diabatic ground state (i.e., in the momentum eigenbasis of the 
uncoupled problem [31]) at t — — 0.5Tb and then propagate it to t — 0.5Tb. The 
survival probability at each time step is obtained by projecting onto the diabatic first 
excited state [3137]. Then the survival probability is averaged over the last 10% of the 
total integration time to estimate the asymptotic value of the survival probability. In 
the end, we average over 100 noise realisations and plot the standard deviation of the 
average survival probability as error bars. 

The data for the noise-driven Wannier-Stark problem are generated in a similar way, 
by integrating the time-dependent Schrodinger equation for a ID or 3D wavefunction. 
At t = 0, the system is prepared in the ground state of the bichromatic lattice 
potential superposed by a small harmonic trap potential (to keep the size of the 
initial wavefunction finite). At t — 0, the harmonic trap potential is disabled and 
replaced by the static force F. The survival probability in the ground band is 
measured by integration over the momentum distribution using appropriate boundaries 
(see [5j|3lJ|l5j[50] f° r more details). While simulations for single-particle wave function 
were performed in one dimension, realistic results for the Gross-Pitaevskii equation 
can only be achieved using simulations in three dimensions (for the used algorithm, 
see [45]). For realistic values of the nonlinearity parameter, we used the values of the 
setup described in [46]: 87 Rb atoms in an optical lattice (di = 426nm) that are initially 
loaded into a cigar-shaped, quasi-one-dimensional optical trap, with 250Hz radial and 
20Hz longitudinal confinement. 
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Appendix B. Unit conversion 

Table Bl: Summary of the new unit system. The table gives relations between our set 
of dimensionless units, Si-units and a set of common experimental units [5|l29|l30]. Here, 
M is the mass of the atoms in the BEC and k L is the wavelength of the laser light used 
to generate the optical lattice. E rec sets the characteristic energy scale of the system. 





Energy 


Momentum 




1 photon exchange 


P _ h ^ k L 

^ rcc — 2M 


Piec = frk L 






Si-units 


dimensionless 


experimental 


Energy 


Esi 


E - Ksi 


E - £sl. 


Time 


tsi 




tcxp = tsi 


Space 


XSI 


x = 2xsik L 


•Eexp — ^Xsi^L 


Force 


Fsi 


Fn - 


TP _ Psitt 
^cxp - E ■ kr 


Potential 


Vsi 


V — 8E 


V — -¥SL. 
' cxp R 
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